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Abstract 

We establish soliton-like asymptotics for finite energy solutions to the Dirac 
equation coupled to a relativistic particle. Any solution with initial state close to 
the solitary manifold, converges in long time limit to a sum of traveling wave and 
outgoing free wave. The convergence holds in global energy norm. The proof uses 
spectral theory and the symplectic projection onto solitary manifold in the Hilbert 
phase space. 
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1 Introduction 

We prove the long time convergence to the sum of a soliton and dispersive wave for the 
Dirac equation coupled to a relativistic particle. The convergence holds in global energy- 
norm for finite energy solution with initial state close to the solitary manifold. Our main 
motivation is to develop the techniques of Buslaev and Perelman [2, 3] in the context of the 
Dirac equation. The development is not straightforward because of known peculiarities 
of the Dirac equation: nonpositivity of the energy, algebra of the Dirac matrices, etc. We 
expect that the result might be extended to nonlinear relativistic Dirac equation relying 
on an appropriate development of our techniques. 

Let ifj(x) G C 4 be a Dirac spinor field in IR 3 , coupled to a relativistic particle with 
position q and momentum p, governed by 



itft(x, t) = [— ia±di — ia 2 d 2 — ia^dz + Pm]tp(x, t) + p{x — q(t)) 
q{t)=p(t)/y/l+j?V)> p(t) = Re(i>{x,t),Vp(x-q(t))) 



xeR* (l.l) 



where p e C(M 3 ,C 4 ) and (•, •) stands for the Hermitian scalar product in L 2 (IR 3 ) ® C 4 . 
Here dj = d/dxj, Oij and /3 are 4x4 Dirac matrices. The standard representation for the 
Dirac matrices aj and j3 (in 2 x 2 blocks) is 



^ a <>={ -I,)' a '=U 0' )• 1 = h2 - 3 (L2) 



where I 2 denotes the unit 2x2 matrix and 



1 \ / -% \ ( 1 



rr ■ i o )■ 01 \ i o )■ a:> > v o -i 

The matrices aj, j — 0, 1, 2, 3 are Hermitian and satisfy the anticommutation relations 



aj, ajQik + a k aj = 2S jk (1.3) 



We will use the following real orthogonality relations 

Pip ■ ajip = 0, j = 1, 3, and a 2 ip ■ ip = 0, ^6l 4 (1.4) 

The system (1.1) is translation-invariant and admits soliton solutions 

s a ,v(t) = (ip v (x - vt - a),vt + a,p v ), p v = v/Vl - v 2 (1.5) 

for all a,i)6K 3 with \v\ < 1. The states S a , v '■= s ajV (0) form the solitary manifold 

S :={S a>v :a,veR 3 ,\v\ < 1} (1.6) 

Our main result is the soliton-type asymptotics 

ip(x, t) ~ ip v± (x - v±t - a±) + W (t)<f>±, t -> ±oo (1.7) 

for solutions to (1.1) with initial data close to the solitary manifold S. Here W (t) 
is the dynamical group of the free Dirac equation, (j>± are the corresponding asymptotic 
scattering states, and the asymptotics hold in the global norm of the Hilbert space L 2 ' 
C 4 . For the particle trajectory we prove that 



q(t) — > v±, q(t) ~ v±t + a±, t — > ±oo 



(1.8) 
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The results are established under the following conditions on the complex valued charge 
distribution: for some v > 5/2 

(1 + \x\Y\d a p\ e L 2 {R 3 ), \a\<3 (1.9) 

We assume p(— x) = p(x), x e R 3 , for the simplicity of calculations. Finally, we assume 
the Wiener condition for the Fourier transform p = (2rr)~ 3 ^ 2 J e lkx p{x)dx 

B(k) =m/3p(k) ■ p(k) > 0, keR 3 (1.10) 

which is the nonlinear version of the Fermi Golden Rule in our case (cf. [4, 13, 14, 15]): the 
nonlinear perturbation is not orthogonal to the eigenfunctions of the continuous spectrum 
of the linear part. The examples are easily constructed. Namely, let us rewrite (1.10) in 
the form 

B(k) = mi\ Pl (k)\ 2 + \p 2 (k)\ 2 - \p 3 (k)\ 2 - |p 4 (A;)| 2 ] > 0, keR 3 (1.11) 

Therefore, we can take e.g. p\ constructed in [12], and p 2 = p% = p^ = 0. 

The system (1.1) describes the charged particle interacting with its "own" Dirac field. 
The asymptotics (1.7)-(1.8) mean asymptotic stability of uniform motion, i.e. "the law 
of inertia". The stability is caused by "radiative damping", i.e. radiation of energy 
to infinity appearing analytically as a local energy decay for solutions to the linearized 
equation. The radiative damping was suggested first by M.Abraham in 1905 in the context 
of the Maxwell-Lorentz equations, [1]. 

One could also expect asymptotics (1.7) for small perturbations of the solitons for 
the relativistic nonlinear Dirac equations and for the coupled nonlinear Maxwell-Dirac 
equations whose solitons were constructed in [6]. Our result models this situation though 
the relativistic case is still open problem. 

Asymptotics of type (1.7)-(1.8) were obtained previously for the Klein-Gordon and 
Schrodinger equations coupled to the particle [8, 10]. More weak asymptotics of type 
(1.7) in the local energy norms, and without the dispersive wave, were obtained in [7] and 
[11] for the Maxwell-Lorentz and wave equations respectively. 

Let us comment on our approach. For ID translation invariant Schrodinger equation, 
asymptotics of type (1.7) were proved for the first time by Buslaev and Perelman [2, 3, 4], 
and extended by Cuccagna [5] for higher dimensions. Here we develop the approach [8] 
where the general Buslaev and Perelman strategy has been developed for the case of the 
Klein-Gordon equation: i) symplectic orthogonal decomposition of the dynamics near 
the solitary manifold, ii) modulation equations for the symplectic projection onto the 
manifold, and iii) the time decay in the transversal directions, etc (see more details in 
Introduction [8]). We prove the asymptotics (1.7)- (1.8) in Sections 3-11 developing this 
general strategy. One of difficulties is caused by well known nonpositivity of the Hamil- 
tonian for the Dirac equation. Respectively, the energy conservation does not provide a 
priori estimate for the solution. We obtain linear in time estimate for L 2 norm of the 
solution using unitarity of the free Dirac propagator. The main novelty in our case is 
thorough establishing the appropriate decay of the linearized dynamics in Sections 12-17, 
and Appendices A, B, and C: 

I. Main difficulty lies in the proof of the decay ~ t~ 3 l 2 in weighted norms for the free 
Dirac equation. Here we prove the decay for the first time (Lemma 17.1). The proof relies 
on the "soft version" of the strong Huygens principle for the Dirac equation. Namely, the 
free Dirac propagator is concentrated mainly near the light cone, while the contribution 
of the inner zone is a Hilbert-Schmidt operator. 
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II. Next difficulty lies in the computation of the spectral properties of the linearized 
equation at the soliton. We do not postulate any spectral properties of the equation in 
contrast to majority of the works in the field. Namely, we find that under the Wiener 
condition (1.10), the discrete spectrum consists only from zero point with algebraic mul- 
tiplicity 6 (Lemma 16.2). The multiplicity is totally due to the translation invariance of 
the system (1.1). 

III. Moreover, we exactly calculate the symplectic orthogonality conditions (16.7) for 
initial data of the linearized equation. These conditions are necessary for the proof of the 
decay. 

IV. All computations differ significantly from the case of the Klein-Gordon equation [8] 
because of the algebra of the Dirac matrices. An important role play the real orthogonality 
relations (1.4) for the Dirac matrices. 

Our paper is organized as follows. In Section 2, we formulate the main result. In Sec- 
tion 3, we introduce the symplectic projection onto the solitary manifold. The linearized 
equation is considered in Sections 4 and 5. In Section 6, we split the dynamics in two 
components: along the solitary manifold, and in transversal directions. The time decay 
of the transversal component is established in sections 7-10 using the time decay of the 
linearized dynamics. In Section 11 we prove the main result. In Sections 12 - 16 we justify 
the time decay of the linearized dynamics relying on the weighted decay for the free Dirac 
equation in a moving frame which is proved in Section 17. In Appendices A, B and C we 
collect some technical calculations. 



2 Main results 

2.1 Existence of dynamics 

We consider the Cauchy problem for the system (3.1) which we write as 

Y(t) = F(Y(t)), tER: Y(0) = Y (2.1) 

Here Y(t) = (ip(t), q(t),p(t)), Y = (ip(0), qo,Po), and all derivatives are understood in the 
sense of distributions. To formulate our results precisely, we need some definitions. We 
introduce a suitable phase space for equation (2.1). Let L 2 a , aeR, denote the weighted 
Agmon spaces with the norm \\ip\\ a = W^Wl'I '■= ||(1 + M) a |^| \\l 2 , where L 2 = L 2 (R 3 ). 

Definition 2.1. i) The phase space £ is the Hilbert space L\ © R 3 © R 3 of states Y = 
(ip,q,p) with the finite norm 

\\Y\\e = U\\o + \q\ + \p\ 
ii) £ a is the space L 2 a © R 3 © IR 3 with the finite norm 

\\y\u = \\Y\\ £a = u\u + \q\ + \p\ 

Proposition 2.2. Let (1.9) hold. Then 

(i) For every Y G £ the Cauchy problem (3.1) has a unique solution Y{t) e C(R, £). 

(ii) For every t G R, the map U(t) : Y i— > Y{t) is continuous on £. 
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Proof. Step i) First, let us fix an arbitrary b > and prove (i)-(ii) for Y G £ such that 
HV'ollo < b and \t\ < e = e{b) for some sufficiently small e(b) > 0. Let us rewrite the 
Cauchy problem (2.1) as 

Y(t) = F 1 (Y(t)) + F 2 (Y(t)), teR: Y(0) = Y (2.2) 

where 

Fi : Y H> ((-ajdj - ipm)ip, 0, 0) 

F 2 : Y i— >• {—ip{x — q), p/ a/1 + p 2 , Re J ■ V p{x — g)rfx) 

The Fourier transform provides the existence and uniqueness of solution Yi(t) G C(R, £) 
to the linear problem (2.2) with F 2 = 0. Let U\(t) : lo l— ^ be the corresponding 

strongly continuous group of bounded linear operators on £ . Then (2.2) for Y(t) G C(R, £) 
is equivalent to the integral Duhamel equation 

t 

Y(t) = U 1 (t)Y + J dsU l (t-s)F 2 (Y(s)) (2.3) 
o 

because F 2 (y(-)) G C(R, £) in this case. The latter follows from local Lipschitz continuity 
of the map F 2 in £: for each b > there exist a x = x(5) > such that for all 
Yi = (ipi,qi,Pi),Y 2 = (i^ 2 ,q 2 ,p 2 ) G £ with HV'illo, H^Ho < b, 

\\F2(Y 1 )-F 2 (Y 2 )\\e<x\\Y 1 -Y 2 \\e 

Therefore, by the contraction mapping principle, equation (2.3) has a unique local solution 
Y(-) G C([— e,e],£) with e > depending only on 6. 

S'tep iij Second we derive a priori estimate. Consider ip G C£°(M 3 ) <g) C 4 . Then 

^IMIo = y (V' • "0 + "0 • = y"(# • ~q)-iip- p(x - q))dx < C\\i/j\\ 
Hence, 

||^(*)l|o<ict+||V(0)|| 

Now, the last two equalities (1.1) imply a priori estimates for \p\ and |g|. The a priori 
estimates for general initial data ipo £ -^o follow by approximating initial data by the 
functions from C^(R 3 ) <g> C 4 . 

Step Hi) Properties (i)-(ii) for arbitrary t G R now follow from the same properties for 
small \t\ and from a priori estimate. □ 

2.2 Solitary manifold and main result 

Let us compute the solitons (1.5). The substitution to (1.1) gives the stationary equations 

-ivjdjip v (y) = [-iaijdj + pm]if> v (y) + p(y) 



v = Pvl v 71 + P% = Re y tp v (y) ■ V ' p{y) dy 

Applying Fourier transform to the first equation in (2.4) we obtain 

(—Vjkj + ajkj — [5m)ip v {k) = p(k) 



(2.4) 
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hence 

7 = fojfy + ~ I3m)p(k) = ( Vj kj + gjfcj - (3m)p{k) 

(vjkj + ctjkj — j3m)(vjkj — atjkj + /3m) k 2 + m 2 — (vjkj) 2 

The soliton is given by the formula 

?7 /" e" m|7(y_:r) ii +(s/ " :E)±l p(|/)(i 3 |/ w 

M*) = ^(vA+*A+iPrn) J Hy _ x)ii + iy _ x)±l , P« = 7^ = (2.6) 

It remains to prove that the last equation of (2.4) holds. Indeed, Parseval identity and 
equality (2.5) imply 

RejMy)-d jP (y)dy = Bejikrf v (k) ■ p(k)dk = Rejikj ^h^h^Mt &1 dk = 

since the integrand is pure imaginary function. Hence, the soliton solution (1.5) exists 
and is defined uniquely for any couple (a,v) with \v\ < 1 and a G M 3 . Let us denote by 
V := {v G M 3 : \v\ < 1}. 

Definition 2.3. A soliton state is S(a) := (ip v (x — b), b, v ), where a := (b,v) with b G M 3 
and v G V . 

Obviously, the soliton solution admits the representation S(cr(t)), where 

a(t) = (b(t),v(t)) = (vt + a, v) (2.7) 

Definition 2.4. A solitary manifold is the set S := {S(a) : a G £ := 1R 3 x V}. 

The main result of our paper is the following theorem. 

Theorem 2.5. Let (1.9), and the Wiener condition (1.10) hold. Let v > 5/2 be the 
number from (1.9), and Y(t) be the solution to the Cauchy problem (2.1) with the initial 
state Y which is sufficiently close to the solitary manifold: 

d :=dist£,(y ,«S) < 1 (2.8) 
Then the asymptotics hold for t — > ±oo, 

q(t) = v± + 0(\t\- 2 ), q(t) = v ± t + a± + 0{\t\- 1 ) (2.9) 

ip(x,t) = ip v ±(x - v±t - a±) + W (t)<t>± + r±{x,t) (2.10) 

with 

\\r ± (t)\\ Q = 0(\t\-^) (2.11) 

It suffices to prove the asymptotics (2.9), (2.10) for t — > +oo since the system (1.1) is 
time reversible. 
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3 Symplectic projection 



3.1 Hamiltonian structure 

Denote ip\ = Retp, ip 2 = Im^, p\ = Rep, p2 = Imp, a 2 = —ia 2 . Then the system (1.1) 
reads 

V>i(x, t) = -(audi + a 3 d 3 )ij 1 (x, t) + (a 2 d 2 + f3m)^ 2 (x, t) + p 2 (x - q(t)) 
ip 2 {x,t) = -(a 2 d 2 + pm)ijji(x,t) - (aidi + a 3 d 3 )ip 2 (x,t) - pi(x - q(t)) 



q(t)=P(t)/V 1 +P 2 ( t ) 
p(t) = J (4> 1 (x,t)-Vp 1 (x-q(t))+ip 2 (x,t)-Vp 2 (x-q(t)))dx 

This is a Hamilton system with the Hamilton functional 



x ER 3 (3.1) 



H(il>i, i>2,q,p) = ^ y (^i • («2^2 + /3^)^i +^2 • (a 2 d 2 + f3m)ip 2 + • (ai9i + a 3 9 3 )^ 2 )^ 
+ J (ip!(x) ■ pi(x - q) + ip 2 (x) ■ p 2 (x - q))dx + ^/T^j? (3.2) 



Equation (3.1) can be written as a Hamilton system 



Y = JDH(Y), Y = {^ 2 ,q,p), J:- 



( J 4 \ 

-J 4 

/ 3 

V -j 3 j 



(3.3) 



where _D"H is the Frechet derivative with respect to ipik, ip 2 k, k — 1,2, 3, 4, p and g of the 
Hamilton functional. 



3.2 Symplectic projection onto solitary manifold 

Let us identify the tangent space to £, at every point, with £. Consider the symplectic 
form (1 defined on £ by Si = / A AfeM dx + dqA dp, i.e. 

si(y 1 ,r 2 ) = (y 1 ,jy 2 ), = e£, j = 1,2 (3.4) 

where 

(Y\ Y 2 } := , + . </> 2 } + q 1 ■ ? 2 + P 1 • V 2 

and = stands for the scalar product or its different extensions. 

It is clear that the form Q is non-degenerate, i.e. 

tt(Y\ Y 2 ) = for every Y 2 E £ ==► F 1 = 

Definition 3.1. i) Y 1 \ Y 2 means that Y 1 E £, Y 2 E £, and Y 1 is symplectic orthogonal 
to Y 2 , i.e. n{Y\Y 2 ) = 0. 

ii) A projection operator P : £ — > £ is called symplectic orthogonal if Y 1 \ Y 2 for 
Y 1 E KerP and Y 2 E ImP. 
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Let us consider the tangent space Ts( a )S to the manifold S at a point S(a). The 
vectors Tj := d aj S(<j), where d aj := and d aj+3 := d Vj with j = 1,2,3, form a basis in 
T a S. In detail, 

Tj = Tj (v) := d b .S(cr) = ( -djijj vl (y) , -dji/j v2 (y) , ej , ) . 12g 
r i+3 = T j+3 (v) := d Vj S(a) = ( d Vj ip vl (y), d Vj ip v2 (y) , 0, <9^p„) 

(3.5) 

where = Re^, ■0„2 = Im^, y :— x— b is the "moving frame coordinate", e\ = (1, 0, 0) 
etc. Let us stress that the functions Tj will be considered always as the functions of y, 
not of x. Formula (2.6) and condition (1.9) imply that 

Tj(v)e£ a , veV, j = l,...,6, VaeR (3.6) 

Lemma 3.2. The matrix with the elements fl(ri(v),Tj(v)) is non- degenerate Vi> G V. 

Proof. The elements are computed in Appendix A. As the result, the matrix Q(ri,Tj) has 
the form 

n(t,) := i<>iv7.v :;!)/ , , „ = ( \ +{v) ^ +{V) ) (3.7) 

where the 3 x 3-matrix Q + (v ) equals 

Q+( v ) = R+(1- v 2 )~ 1/2 E + (1 - v 2 )~ 3/2 v ® v (3.8) 
Here K is a symmetric 3 x 3-matrix with the elements 

where i3(/c) > is defined in (1.10). The matrix K is the integral of the symmetric 
nonnegative definite matrix k <8> k = (kikj) with a positive weight. Hence, the matrix K 
is nonnegative definite. Since the unite matrix E is positive definite, the matrix Q + (v) is 
symmetric and positive definite, hence non-degenerate. Then the matrix Q(ti,tj) also is 
non-degenerate. □ 

Let us introduce the translations T a : (ip(-),q,p) i-> (*0(- — a),g + a,p), a G M 3 . 
Note that the manifold S is invariant with respect to the translations. Let us denote 
v (p) ■= p/y/l +p 2 for p G 1R 3 . 

Definition 3.3. i) For any a G R and v < 1 denote by £ a (v) = {Y = (vp,q,p) G : 
|u(p)| < v}. We set £(v) := S (v). 

ii) For any v < 1 denote by £({;) = {a = (6, v ) : b G 1R 3 , |f| < -u}. 

The next Lemma provide that in a small neighborhood of the soliton manifold S a 
"symplectic orthogonal projection" onto S is well-defined. 

Lemma 3.4. (cf.[8, Lemma 3.4]) Let (1.9) hold, a G R andv < 1. Then 

i) there exists a neighborhood O a (S) of S in S a and a map II : O a (S) — > S such that II 
is uniformly continuous on O a (S) n£ a (v) in the metric of £ a , 

UY = Y for YeS, and Y-S\T S S, where S = UY (3.10) 

ii) O a (S) is invariant with respect to the translations T a , and 

UT a Y = T a UY, for Y G O a {S) and a G R 3 

Hi) For any v < 1 there exists a v < 1 s.t. IIY = S(a) with a G S({i) for Y G 
O a (S)D£ a (v). 

iv) For any v < 1 there exists an r a {y) > s.t. S(a) + Z G O a (S) if a G S(t>) and 
\\Z\\ a < r a (v). 
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We will call II a symplectic orthogonal projection onto S. 

Corollary 3.5. The condition (2.8) implies that Y = S + Z where S = S(a ) = UY , 
and 

||^o||,<l (3.11) 

Proof. Lemma 3.4 implies that II Yo — S is well defined for small do > 0. Furthermore, the 
condition (2.8) means that there exists a point Si E S such that \\Y — Si\\ u = d . Hence, 
Y , Si E O u (S) P\£ u {v) with a v < 1 which does not depend on d for sufficiently small do- 
On the other hand, II Si = Si, hence the uniform continuity of the map II implies that 
||Si — S\\ v — > as d — > 0. Therefore, finally, ||Zo||„ = ll^o - S\\ u < \\Y — Si\\ u + \\Si — S\\ u < 
d + o(l) 1 for small do- □ 



4 Linearization on solitary manifold 

Let us consider a solution to the system (3.1), and split it as the sum 

Y(t) = S(a(t)) + Z(t) 



(4.1) 



where a(t) = (b(t),v(t)) G E is an arbitrary smooth function of t G R. In detail, denote 
Y = (ip,q,p) and Z = (^,g,P). Then (4.1) means that 

ii(x,t)=iP v{t) (x-b(t)) + V(x-b(t),t), q(t)=b(t)+Q(t), p(t)=p v{t) + P(t) (4.2) 

Let us substitute (4.2) to (1.1), and linearize the equations in Z. Setting y = x — b(t) 
which is the "moving frame coordinate" , we obtain that 

i> = v-V v ip v (y)-b-ViP v (y) + i>(y : t)-b-V'S>(y,t) 
= l-ajdj - tf3m]{^ v (y) + *(y,t)) - ip{y - Q) 

q = b + Q = - 

y/1 + ( Pv + py 

P = v-V v p v + P = Re(i> v (y) + y(y,t),Vp(y-Q)) 

Let us extract linear terms in Q. First note that p(y — Q) = p(y) — Q ■ V p(y) + Ni(Q), 
Vp{y-Q)=V P {y)-V{Q-V P {y)) + Ni{Q). 

The condition (1.9) implies that for Ni(Q) and Ni(Q) the bound holds, 



(4.3) 



1^)11, + 11^(3)11, <a(Q)Q 2 



(4.4) 



uniformly in \Q\ < Q for any fixed Q, where v is the parameter from Theorem 2.5. Second, 
the Taylor expansion gives 

h = = m =v + kP~ v(v ■ P)) + N 2 (v, P) 
VI + (p v + P) 2 7 



where 1/7 = Vl — v 2 = (1 + pi) and 

\N 2 (v,P)\ < C{v)P 2 



(4.5) 
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uniformly with respect to \v\ < v < 1. Using the equations (2.4), we obtain from (4.3) 
the following equations for the components of the vector Z(t): 



(4.6) 



%, f) = [-ajdj - ipmYi!(y, t) + b ■ V*(y, t) + iQ- Vp{y) 

+ (b-v)- vMv) - * ■ VvMv) - 

Q(t) = ±{E-v® v)P + (v-b) + N 2 

p(t) = -v. V vPv + Re{V(y,t),Vp(y)) + Re{VMv),Q -Vp(y)) + N 3 (v, Z) 

where N 3 (v, Z) = - Re(V^, JVi(Q))-Re<# , V(Q-Vp))+Re(^, Ni(Q)). Clearly, N 3 (v, Z) 
satisfies the following estimate 



\N 3 (v,Z)\<C v (p,v,Q) Q 2 +\\V\\-M 



(4.7) 



uniformly in \v\ < v and \Q\ < Q for any fixed v < 1. For the vector version Z 
(*i, * 2 , Q, P) with = Re ^ 2 = Im * we rewrite the equations (4.6) as 



Z{t) = A{t)Z{t)+T{t) + N{t), te 



(4.8) 



Here the operator A(t) = A v>w (t) depends on two parameters, v = v(t), and w = b(t) and 
can be written in the form 







^2 




i Q 









( —Q.\d\— a 2 d 2 +w V a 2 d 2 + /3m 

-(af 2 <9 2 + /3m) -«i<9i-Q;2i92+w -V 


(-,v Pl ) 



-Vp 2 " 





<-,Vp2 



\ 













^2 




By 


















(4.9) 


a, Z) in (4.8) 



stand for 



J- 7 



/ (u> - v ) • V^i - v ■ V v 1p v i \ 

(w - V ) • V^2 - • V„V'r,2 
V — W 

-v ■ V v p v 



N(a,Z) 



( N l2 {Z) \ 
-N n (Z) 
N 2 (v,Z) 
V N 3 (v,Z) J 



where v = v(t), w = w(t), a = a(t) = (b(t),v(t)), and Z 
(4.5) and (4.7) imply that 



(4.10) 

Z(t). The estimates (4.4), 
(4.11) 



\\N(a,Z)\\ v <C(v,Q)\\Z\\ 2 _ v 
uniformly in a E £(£>) and || ^|| _^ < r_ v (v) for any fixed v < 1. 

Remark 4.1. i) The term in the right hand side of the equation (4.8) is linear 

in Z(t), and N(t) is a high order term in On the other hand, T(t) is a zero order 

term which does not vanish at Z(t) = since S(a(t)) generally is not a soliton solution if 
(2.7) does not hold (though S(a(t)) belongs to the solitary manifold), 
ii) Formulas (3.5) and (4.10) imply: 



T(t) = -^2[(w- v)iT t + VlT l+3 \ 



(4.12) 



and hence T(t) E Ts(a(t))<S, t E 
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5 Linearized equation 



Here we collect some Hamiltonian and spectral properties of the generator (4.9) of the 
linearized equation. First, let us consider the linear equation 



X(t) = A VjW X(t), t el, v e V, we 



(5.1) 



Lemma 5.1. (cf. Lemma 5.1 [8]) i) For any v G V and w G M. 3 the equation (5.1) can 
be written as the Hamilton system (cf. (3.3)), 



X(t) = JDH v , w (X(t)), t G 



(5.2) 



where D% VyW is the Frechet derivative with respect to ^!ik, ^2k, k = 1,2,3,4, P and Q of 
the Hamilton functional 

U V , W {X) = X - J(^ 1 ■ (a 2 d 2 + (3m)m 1 + ^ 2 • (a 2 d 2 + (3m)m 2 + 2^ • (a^ + a 3 d 3 )V 2 )dy 

+ J pAv)Q-V*idy + \p-B v P-±(Q-V^ v M,Q-VpM), X = (* 1 ,* 2 ,Q,P)eS 

(5.3) 

ii) The skew-symmetry relation holds, 

n(A V)W x u x 2 ) = -q(x 1 ,a VjW x 2 ), x l es, x 2 e h\r 3 ) ®h\r 3 ) ©m 3 ©r 3 (5.4) 

Lemma 5.2. The operator A VtW acts on the tangent vectors Tj(v) to the solitary manifold 
as follows, 

A VtW [T j (v)} = (w - v) ■ Vtj(u), A VjW [r j+3 (v)] = (w - v) ■ Vr i+3 (w) + Tj(v), j = 1, 2, 3 

(5.5) 



Proof. In detail, we have to show that 

( -djip vl \ 
-djip v2 



/ (v-w)- Vdjip vl \ 
v -w) ■ Vdjip v2 




A r 



\ 







) 


V o 


( d Vj ijj vl \ 




( 


(w 


- v) ■ Vd Vj ip vl \ 


d Vj ip v2 






(w 


- v) ■ Vd Vj ip v2 














V d VjPv J 











/ -djijj vl \ 
V o J 



(5.6) 



Indeed, differentiate the equations (2.4) in bj and Vj, and obtain that the derivatives of 
soliton state in parameters satisfy the following equations, 

— v ■ Vdjipv = [~ ol • V — if3m]djip v — idjp 
-djip v - v ■ S7d Vj ip v = [-a ■ V - iPm]d Vj il> v 



d Vi p v = e,{l-v 2 Y l l 2 + 



f A _ v 2y/2 

= (d Vj ip vl , Vpi) + (d Vj ip v2 , Vp 2 ) 
for j = 1, 2, 3. Then (5.6) follows from (5.7) by definition of A in (4.9) 



(5.7) 



□ 
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Corollary 5.3. Let w = v e V. Then Tj(v) are eigenvectors, and Tj +3 (v) are root vectors 
of the operator A v>v , corresponding to zero eigenvalue, i.e. 

Av,v[Tj( v )\ = °. A vA t j+3( v )} = Tj(v), j = 1,2,3. (5.8) 

Remark 5.4. For a soliton solution of the system(3.1) we have b — v, v — 0, and hence 
T(t) = 0. Thus, the equation (5.1) is the linearization of the system (3.1) on a soliton 
solution. In fact, we do not linearize (3.1) on a soliton solution, but on a trajectory 
S(a(t)) with a(t) being nonlinear in t. We will show later that T(t) is quadratic in Z(t) 
if we choose S(a(t)) to be the symplectic orthogonal projection of Y(t). Then (5.1) is 
again the linearization of (3.1). 

6 Symplectic decomposition of dynamics 

Here we decompose the dynamics in two components: along the manifold S and in 
transversal directions. The equation (4.8) is obtained without any assumption on a(t) in 
(4.1). We are going to choose S(a(t)) := UY(t), but then we need to know that 

Y(t) e 0-„(S), teR (6.1) 

It is true for t — by our main assumption (2.8) with sufficiently small d > 0. Then 
S{a(0)) = ny(0) and Z(0) = Y(0) - S(a(0)) are well defined. We will prove below that 
(6.1) holds if d is sufficiently small. Let us choose an arbitrary v such that \v (0)| < v < 1 
and let 8 — v — \v(0)\. Denote by r^ v (v) the positive numbers from Lemma 3.4 iv) 
which corresponds to a = —v. Then S(a) + Z e 0^ U (S) if a — (b,v) with \v\ < v 
and < r- v (v). Note that ||Z(0)||_ V < r_ v {v) if d is sufficiently small. Therefore, 

S(a(t)) = UY(t) and Z(t) = Y(t) - S(a(t)) are well defined for t > so small that 
\v | < v and ||^(/:) < r- v (v). This is formalized by the following standard definition. 

Definition 6.1. t* is the "exit time", 

U = sup{* > : 11^(5)11-^ < r_„(u), \v(s) - v(0)\ < 5, < s < t} (6.2) 

One of our main goals is to prove that t* = oo if d is sufficiently small. This would 
follow if we show that 

\\Z(t)\\- v <r- v (v)/2, \v(s)-v(0)\<6/2, 0<t<U (6.3) 

Note that 

\Q(t)\ <Q:= r- v (v), 0<t<U (6.4) 
Now N(t) in (4.8) satisfies, by (4.11), the following estimate, 

\\N(t)\\u <C v (v)\\Z(t)\\i v , 0<t<U (6.5) 

6.1 Longitudinal dynamics: modulation equations 

From now on we fix the decomposition Y(t) = S(a(t)) + Z(t) for < t < by setting 
S(a(t)) = HY(t) which is equivalent to the symplectic orthogonality condition of type 
(3.10), 

Z(t)\T S (a(t))S, 0<t<U (6.6) 

This allows us to simplify drastically the asymptotic analysis of the dynamical equations 
(4.8) for the transversal component Z(t). As the first step, we derive the longitudinal 
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dynamics, i.e. the "modulation equations" for the parameters cr(t). Let us derive a 
system of ordinary differential equations for the vector <r(t). For this purpose, let us write 
(6.6) in the form 

n(Z(t),Tj(t)) = 0, j = 1,...,6, 0<t<U (6.7) 

where the vectors Tj(t) = Tj(a(t)) span the tangent space Ts( a {t))S- Note that a(t) = 
(b(t),v(t)), where 

\v(t)\ < v < 1, 0<t<U (6.8) 

by Lemma 3.4 iii). It would be convenient for us to use some other parameters (c,v) 
instead of a — (b,v), where 

c(t) = b(t) - [ v(r)dr, c(t) = b(t) - v(t) = w(t) - v(t), < t < U (6.9) 
Jo 

The following statement can be proved similar to the Lemma 6.2 in [8]. 

Lemma 6.2. Let Y(t) be a solution to the Cauchy problem (3.1), and (4-1), (6-7) hold. 
Then 

\c(t)\ + \v(t)\<C(v)\\Zt u (6.10) 
6.2 Decay for transversal dynamics 

In Section 11 we will show that our main Theorem 2.5 can be derived from the following 
time decay of the transversal component Z(t): 

Proposition 6.3. Let all conditions of Theorem 2.5 hold. Then t* = oo, and 

We will derive (6.11) in Sections 7-10 from our equation (4.8) for the transversal 
component Z(t). This equation can be specified using Lemma 6.2. Indeed, the lemma 
implies that 

\\T(t)\\ v <C(v)\\Z(t)\Zv 0<t<U (6.12) 

by (4.10) since w — v — c. Thus (4.8) becomes the equation 

Z(t) = A(t)Z(t) + N(t), 0<t<U (6.13) 

where A(t) = A v ( t ),w(t), and N(t) := T(t) + N(t) satisfies the estimate 

\\N(t)\\, < C{v,Q)\\Z{t)\\ 2 -v 0<t<U (6.14) 

In all remaining part of our paper we will analyze mainly the basic equation (6.13) to 
establish the decay (6.11). We are going to derive the decay using the bound (6.14) and 
the orthogonality condition (6.6). 

Similarly [8] we reduce the problem to the analysis of the frozen linear equation, 

X(t)=A 1 X(t), teR (6.15) 

where A\ = A VltVl with v\ = v(ti) and a fixed t\ e [0,t*). Then we can apply some 
methods of scattering theory and then estimate the error by the method of majorants. 
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Note, that even for the frozen equation (6.15), the decay of type (6.11) for all solutions 
does not hold without the orthogonality condition of type (6.6). Namely, by (5.8) the 
equation (6.15) admits the secular solutions 

3 3 

x(t) = c jtj(v) + D A T A v ) t + ( 6 - 16 ) 

1 1 

which arise by differentiation of the soliton (1.5) in the parameters a and v in the moving 
coordinate y = x — v\t. Hence, we have to take into account the orthogonality condition 
(6.6) in order to avoid the secular solutions. For this purpose we will apply the corre- 
sponding symplectic orthogonal projection which kills the "runaway solutions" (6.16). 

Remark 6.4. The solution (6.16) lies in the tangent space Ts^S with o\ = (&i,i>i) (for 
an arbitrary bi G K) that suggests an unstable character of the nonlinear dynamics along 
the solitary manifold. 

Definition 6.5. i) For v G V , denote by Yl v the symplectic orthogonal projection of £ 
onto the tangent space Ts(p)S, and P v — I — 

ii) Denote by Z v = P v £ the space symplectic orthogonal to Ts(a)S with a = (b,v). 
Note that by the linearity, 

n v z = n^v^ivMniv), z), ze£ (6.17) 

with some smooth coefficients Ilji(v). Hence, the projector 11^, in the variable y = x — b, 
does not depend on b, and this explains the choice of the subindex in and P„. 
Now we have the symplectic orthogonal decomposition 

£ = T S (a)S + Z v , a=(b,v), (6.18) 

and the symplectic orthogonality (6.6) can be written in the equivalent forms, 

n w(t) Z(t) = 0, P v{t) Z{t) = Z{t), 0<t<U (6.19) 

Remark 6.6. The tangent space Ts(a)<S is invariant under the operator A VjV by Lemma 
5.3 i), hence the space Z v is also invariant by (5.4): A VtV Z G Z v for sufficiently smooth 
Z G Z v . 

Below in section 12-18 we will prove the following proposition which will be one of the 
main ingredients for proving (6.11). Let us consider the Cauchy problem for the equation 
(6.15) with A = A V;V for a fixed v G V. Recall that the parameter v > 5/2 is also fixed. 

Proposition 6.7. Let the conditions (1.9)- (1.10) hold, \v\ < v < 1, and Xq G £. Then 

i) Equation (6.15), with A = A VfV , admits the unique solution e At X := X(t) G C(R, £) 
with the initial condition X(0) = X . 

ii) For X G Z v D £ v , the decay holds, 



(6.20) 



15 



7 Frozen transversal dynamics 

Now let us fix an arbitrary t\ G [0, £*), and rewrite the equation (6.13) in a "frozen form" 
Z(t)=A 1 Z(t) + (A(t)-A 1 )Z(t) + N(t), 0<t<U (7.1) 
where A 1 = A v(tlhv ( tl) and 

\ 




A(t)-A 1 = 



([w(t)- u(ti)] • V 

[iu(t)-' u(*i)] • V 


o o (v(Vyt)j--0«(ti)j), Vpj) 



B v ( t )—B vi (t) 



/ 

The next trick is important since it allows us to kill the "bad terms" [tu(i) — f(£i)] • V in 
the operator A(t) — A x . 

Definition 7.1. Let us change the variables (y,t) i— >■ (yi,t) = (y + di(t),t), where 

di(t) := / (w(s) - u(ti))ds, < t < *! (7.2) 

Next define 

Z x (t) := (Mvi ~ di(t),t), V 2 (Vi ~ di(t),t), Q(t), P(t)) (7.3) 
Then we obtain the final form of the "frozen equation" for the transversal dynamics 



Zi(t) = AiZi(t) + B^Z^t) + iVi(t), < t < ti 
where Aq(i) = iV(£) expressed in terms of y = y\ — <ii(t), and 



(7.4) 



Bi(t) = 






B v ( t ) - B Vl{t ) 



\ 



1 





\0 {VtyvM-ipv^Vpj) 

Let us estimate the "remaining terms" Bi(t)Zi(t) and Ni(t). 
Lemma 7.2. TTie bound holds 

WBiWZ^W,, < C(v)\\Z(t)\\- v £ \\Z(s)\\ 2 _ u ds, 0<t<h 

Proof. Lemma 6.2 implies 

t t 
\B v[t ) ~B Vl{t) \ < | J v(s)-V v B v{s) ds\ < C(v) J \\Z(s)\\ 2 _ v ds 
ti ti 

similarly, 



(7.5) 



\(V^ vW -^ v{tl)j ),Vp 3 \<C(v) j \\Z{s)\\iJs 

ti 

Therefore, 

11^(^)^(011, = KV(^ (t)j -^(tOjO.Vp^QiWI + |(S w(t) - J B, l(t) )P 1 (0| 

t 

< C(i;)(|Q(0| + |P(0l) / 11^)111^ <C7(i;)||Z(0||-, f 



□ 
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Lemma 7.3. The bounds hold 



\\Ni(t)\l < C(v,Q)(l + IdiWDl^WII-*, < t < h (7.6) 
Proof. For any and d e M 3 we have 

= /l<f(y-^)l 2 (i + bl) 2a ^ = /l<f(i/)l 2 (i + b + rf|) 2 ^ 

< J |<%)| 2 (1 + |y|) 2 "(l + \d\) 2a dy < (1 + Ml) 2a ||$|| 2 , «eK 
Hence, the bound (7.6) follows. □ 

8 Integral inequality 

The equation (7.4) can be written in the integral form: 

Z x (t) = e Al %(0) + f e Al{t - s) [BiZi(s) + iVi(s)]ds, < t < h (8.1) 

Now we apply the symplectic orthogonal projection Pi : = P^) to both sides of (8.1): 

PiZi(f) = e^'P^O) + / e A ^P 1 [B 1 Z 1 {s) + N^ds 

Jo 

The projector Pi commutes with the group e Alt since the space Z\ : = P]£ is invariant 
with respect to e Alt by Remark 6.6. Applying (6.20) we obtain that 



i PlZlW i U < E^ +0J r 



|Pi[£iZi(s) + JV 1 (s)]|| I/ ds 



(l + |t-s|) 3 / 2 

The operator P x = I — IIx is continuous in £ v by (6.17). Hence, (7.5)-(7.6) imply 



+ g(di(t))j[' (1 + |t 1 _ 8|)!l/2 [||ZMI|-^"||Z(t)|| 2 .^t+||Z(5)||^ 



where := sup 0<s<t |di(s)|. Let us introduce the "majorant" 



ds, 0<t< h. 

(8.2) 



m(t) := sup (l + s) 3/2 ||Z(s)||_„ te[0,t,). (8.3) 
se[o,t] 

Now we reduce further the exit time. Denote by e < 1 a fixed positive number which we 
will specify below. 

Definition 8.1. t'^ is the exit time 

< = sup{t G [0, U) : m(s) < e, 0<s<t}. (8.4) 
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To estimate di(t), note that 

w(s) -v(ti) = w(s) -v(s) + v(s) -v(h) = c(s) + J v(r)dr (8.5) 
by (6.9). Hence, (7.2), Lemma 6.2 and the definition (8.3) imply that for t\ < t'^ 



\d l {t)\ = \ f (w(s) - v(t!))ds\ < [ (\c(s)\+ [ \v(T)\dr)ds 

Jti Jt V Js 



< C(u)m 2 (ti) 



/ ' (iTTW + 1 ' (TTrp) ds £ c(i)m2(tl) s c(6) ' 5 ' 5 (1 

(8.6) 



Now we can to replace C(d\) with C(t5) in (8.2): for t\ < 



+ C Wj (l + |/_ 8 |)3/2 pWH-/ 1 11^)11-^+11^)11-, 



*i 

' ds, < t < h (8.7) 



9 Symplectic orthogonality 

Finally, we are going to change Pi-Zi(i) by Z(t) in the left hand side of (8.7). We will 
prove that it is possible using again that d <C 1 in (2.8). 

Lemma 9.1. (cf.[8j) For sufficiently small e > 0, we have for t\ < 

\\Z(t)\\- v < CWPiZ^W-,,, 0<t<h, (9.1) 

where C depends only on p and v. 

Proof. Since {d^t)] < C for t < t± < < then ||Z(t)||-i, < C||£i(*)||-i>, and it suffices to 
prove that 

||Zi(*)||_ 1/ <2||PiZi(*)||-„, < t < t x . (9.2) 
Recall that P±Zi(t) = Zi(t) — H v ( tl )Zi(t). Then estimate (9.2) will follow from 

Wn^z^t)]]^ < ^^(Oll-,, o<t<h. (9.3) 

Symplectic orthogonality (6.19) implies 

n vm z l {t) = o, te[o,h], (9.4) 

where Tl v (t),iZi(t) is Il v ^Z(t) expressed in terms of the variable y± = y + d\(t). Hence, 
(9.3) follows from (9.4) if the difference LL^) — U v u\i is small uniformly in t, i.e. 

||n„ (tl) - n„ (t);1 || < 1/2, o < t < h. (9.5) 

It remains to justify (9.5) for small enough e > 0. Formula (6.17) implies 

ii^Mt) = J^njMtyTjMminMt))^^)), (9.6) 
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where Tj t i(v(t)) are the vectors Tj(v(t)) expressed in the variables y\. Since |g?i(£)| < C 
and Vrj are smooth and fast decaying at infinity functions, then 



ho,Mt)) - TMt))\\» < c\di(t)\ u <c, o<t<u 

for all j = 1, 2, . . . , 6. Furthermore, 



(9.7) 



Tj(v(t))-TMtl)) 



v(s) ■ V v Tj(v(s))ds 



and therefore 



Similarly, 



\Tj(v(t)) -Tjivih)) ||„ < cj \v(s)\ds, 0<t<h 



in^M*)) -!!,■,(«(*!)) | = | J 1 v{s)-V v TL jl {v{s))ds\ < C J * \v{s)\ds, < t < h 



(9.8) 



(9.9) 



since |V„nj;(v(s))| is uniformly bounded by (6.8). Hence, the bounds (9.5) will follow 
from (6.17), (9.6) and (9.7)-(9.9) if we establish that the integral in the right hand side 
of (9.8) can be made as small as we please by choosing e > small enough. Indeed, 



J \v(s)\ds <Cm 2 (ti) J 



ds 



< Ce 2 , < t < h 



(9.10) 
□ 



10 Decay of transversal component 



Here we prove Proposition 6.3. 

Step i) We fix < e < 1 and t'^ = t'^(e) for which Lemma 9.1 holds. Then the bound of 
type (8.7) holds with ||PiZi(t)||_, in the left hand side replaced by ||Z(t) : 



\\Z(t)\\-„< 



c 



Jo 



(1+t) 3 / 2 
1 



(l + | t _ s |)3/2 



ll^(0)||, 

\\Z(s)\\_ v £\\Z(r)\\ 2 _ v dT+m 



ds, < t < tx (10.1) 



for ti < t*. This implies an integral inequality for the majorant m(t) defined in (8.3). 
Namely, multiplying both sides of (10.1) by (1 + t) 3 ^ 2 , and taking the supremum in 
t e [0, ti], we get 



m(*i) < C\\Z(0)\\ v +C sup 



(l + tf/ 2 



mis) 



Is (1 



m 2 (r)dr m 2 (s) 



(l + s )3/ 2 y s (i + r )3 (i + s )3 



ds 



tem] Jo (l + \t-s\f/ 2 
for t\ < i*. Taking into account that m(t) is a monotone increasing function, we get 

mfa) < C\\Z{Q)\\ V + C[m 3 (*i) + m 2 ^)]/^), t, < < (10.2) 

where 



r* (i + tf/ 2 



dr 



1 



(l + s) 3 / 2 J s (1 + r) 3 (i + s y 



ds < I < oo 
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Therefore, (10.2) becomes 

m{ti) <C\\Z(0)\\ u + C7[m 3 (t 1 )+m 2 (t 1 )}, h < < (10.3) 
This inequality implies that m(ti) is bounded for t x < t'^, and moreover, 

m(tx) < dllE'CO)^, tx<< (10.4) 

since m(0) = ||Z(0)|| V is sufficiently small by (3.11). 

Step ii) The constant C\ in the estimate (10.4) does not depend on t* and t'^ by Lemma 
9.1. We choose d in (2.8) so small that 11^(0)^ < e/(2C 1 ). It is possible due to (3.11). 
Then the estimate (10.4) implies that t' # = i* and therefore (10.4) holds for all ti < £*. 
Further, 

t t 
\v(t)-v(0)\< J \v(s)\ds<Cm 2 (t) J < Cm 2 it) 



Hence the both inequalities (6.3) also holds if \\Z(0) \\ u is sufficiently small by (8.3). Finally, 
this implies that t* = oo, hence also t'^ = oo and (10.4) holds for all t\ > if d is small 
enough. It complete the proof of Proposition 6.3. 



11 Soliton asymptotics 

Here we prove our main Theorem 2.5 under the assumption that the decay (6.11) holds. 
First we will prove the asymptotics (2.9) for the vector components, and afterwards the 
asymptotics (2.10) for the fields. 

Asymptotics for the vector components. From (4.3) we have q = b + Q, and from 
(6.13), (6.14), (4.9) it follows that Q = P + 0(\\Z\\ 2 _ V ). Thus, 

q = b + Q = v(t) + c{t) + P{t) + 0{\\Z\\ 2 _ V ) (11.1) 

Bounds (6.10) and (6.11) imply that 

W)\ + \m < ^^yr 1 . * > o (ii.2) 

Therefore, c(f) = c+ + 0(t~ 2 ) and v(t) = v + + C(t" 2 ), t ->■ oo. Since \P\ < \\Z\\_ V , the 
estimate (6.11), and (11.1)-(11.2), imply that 

q(t) = v + + 0(rV 2 ) (11.3) 

Similarly, 

b(t) = c(t)+ I v{s)ds = v + t + a + + Oit- 1 ) (11.4) 
Jo 

hence the second part of (1.8) follows: 

q(t) = b(t) + Q(t) = v + t + a + + C(r x ) (11.5) 
since Q(t) = C(t~ 3 / 2 ) by (6.11). 

Asymptotics for the fields. For the field part of the solution i/)(x,t) let us define the 
accompanying soliton field as ip v (t)(x — q(i)), where we define now v(t) = q(t), cf. (11.1). 
Then for the difference z(x,t) = ip(x,t) — ip v ( t )(x — q(t)) we obtain the equation 

z(x, t) = [—eejdj — ij3m]z(x, t) — iv ■ V v ^ v (t) ( x ~ q(t)) 
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Then 

z(t) = W o (t)z(0) - [ W (t - s)\iv(s) • VWv (s) (- - q(s))]ds (11.6) 
Jo 

To obtain the asymptotics (2.10) it suffices to prove that z(t) = Wo(t)(f> + + r+(t) with 
some 4> + G L\ and ||r + (t)|| = Oiir 1 ! 2 ). This is equivalent to 

W (-t)z(t) = <f> + + r' + (t) (11.7) 

where ||r+(t)||o = 0{ir l l 2 ) since VFo(i) is a unitary group in Lq by the charge conservation 
for the free Dirac equation. Finally, (11.7) holds since (11.6) implies 

W (-t)z(t) = z(0) - [ W (-s)f(s)ds, f(s) = iv(s) • V v ^ v(i) (- - q(s)) 

Jo 

where the integral in the right hand side converges in Lq with the rate 0(t~ l l 2 ). The 
latter holds since 1 1 VKo ( — s)/(s)||o = 0(s~ 3 / 2 ) by the unitarity of Wo(—s) and the decay 
rate ||/(s)||o = 0(s~ 3 / 2 ). Let us prove this rate of decay. It suffices to prove that 
\v(s)\ = 0(s- 3/2 ), or equivalents \p(s)\ = C(s" 3/2 ). Substitute (4.2) to the last equation 
of (1.1) and obtain 

p(t) = Re J [ij v(t) (x - b(t)) + *(x - b(t),t)] Vp{x - b(t) - Q(t))dx 

= Re / i; v{t) (y)Vp(y)dy + Re / ^ v(t) (y) [V p{y - Q(t)) - Vp(y)\ dy 

+ 



Re J *{y,t)Vp{y-Q{t))dy 



The first integral in the right hand side is zero by the stationary equations (2.4). The 
second integral is (D(t~ 3 / 2 ), since Q(t) = (D(t~ 3 / 2 ), and by the conditions (1.9) on p. 
Finally, the third integral is 0(t~ 3 ^ 2 ) by the estimate (6.11). The proof is complete. 



12 Decay for linearized dynamics 

In remaining sections we prove Proposition 6.7. Here we discuss our general strategy of 
the proof. We apply the Fourier-Laplace transform 

POO 

X(\)= e- xt X(t)dt, Re A > (12.1) 
Jo 

to (6.15). According to Proposition 6.7, we expect that the solution X(t) is bounded in 
the norm || • ||_„. Then the integral (12.1) converges and is analytic for ReA > 0. We 
will write A and v instead of A\ and v\ in all remaining part of the paper. After the 
Fourier-Laplace transform (6.15) becomes 

AX(A) = AX"(A) + X , Re A > (12.2) 

Let us stress that (12.2) is equivalent to the Cauchy problem for the functions X{t) G 
Cb([0, oo); £- v ). Hence the solution X{t) is given by 

X(\) = -(A - A)" 1 X , Re A > (12.3) 

if the resolvent R(X) = (A — A)" 1 exists for ReA > 0. 
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Let us comment on our following strategy in proving the decay (6.11). First, we will 
construct the resolvent R(X) for Re A > and prove that it is a continuous operator in 
£_ v . Then X(X) G S- u and is an analytic function for Re A > 0. Second, we have to 
justify that there exist a (unique) function X(t) G C([0, oo); £~ u ) satisfying (12.1). 

The analyticity of X(X) and Paley- Wiener arguments (see [9]) should provide the 
existence of a £_ v - valued distribution X(t), t G R, with a support in [0, oo). Formally, 

A- X X = X(t) = — [ e iwt X{iu + 0)du, t G R (12.4) 
2tt Jr 

However, to check the continuity of X(t) for t > 0, we need additionally a bound for 
X{iuj + 0) at large \oj\. Finally, for the time decay of X(t), we need an additional infor- 
mation on the smoothness and decay of X(ko + 0). More precisely, we should prove that 
the function X(iu + 0) 

i) is smooth outside u> = and oj = ±/x, where \i = /i(v) > 0, 

ii) decays in a certain sense as \u\ — > oo. 

iii) admits the Puiseux expansion at oj — ±fi. 

iv) is analytic at oj — if X G Z v := P v £ and X G £ u . 

Then the decay (6.11) would follow from the Fourier-Laplace representation (12.4). 



13 Solving the linearized equation 



Here we construct the resolvent as a bounded operator in £_ v for Re A > 0. We will write 
(*i,*2,<5,-P) instead of (*i(y, A), ^ 2 (y, A), Q(X), P(A)) to simplify the notations. Then 
(12.2) reads 



(A-X) 



# 2 

Q 

\ P ) 



#02 

Qo 

V P o ) 



It is the system of equations 

- a 3 d 3 + v ■ V - A)*i + (/3m + a 2 <9 2 )\i> 2 - Q • Vp 2 = -*oi 



-(/3m + a 2 9 2 )^i + (-ai9i - a 3 <9 3 + u • V - A)^ 2 + Q ■ Vpi = 



02 



P„P - AQ = -g 



-(Vfj, pj) + (V^-, q • Vp,) - AP = -P 



13.1) 



Sfep Let us study the first two equations. First, we compute the matrix integral kernel 
G\(y — y') of the Green operator 



G x 

In Fourier space 
G x (k) = 



—a\d\ — a 3 <9 3 + v ■ V — A 
-/3m - a 2 d 2 



(5m + a 2 d 2 
—a\d\ — a 3 <9 3 + v ■ V — A 



(13.2) 



iotiki + ict^kz — ivk — A (3m — a 2 k 2 

—(5m + a 2 k 2 iot\k\ + ia 3 k 3 — ivk — X 



, v k = ^v j k j 

3=1 
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To invert the matrix, let us solve the system 

(13.3) 



afi + 6/2 = 91 
-bfi + a/ 2 = 02 



where a = iai/ci + ^0:3^3 — iv k — A, 6 = (3m — a 2 k 2 . Multiplying the first equation of (13.3) 
by c = —iaiki — ia 3 k 3 — ivk — A and the second equation by —b, we obtain 



ca/i + cbf 2 = cg x 
b 2 fi - cbf 2 = -bg 2 



(13.4) 



since ba = cb by the anticommutations (1.3). Further, b 2 + ac = k 2 + m 2 + (u>&; + A) 2 . 
Therefore, summing up the equations (13.4), we obtain that 

, cgi - bg 2 

Ji 



k 2 + m 2 + (ivk + A) 2 
Similarly, we obtain 

651 + cg 2 

h 



k 2 + m 2 + (?t>/c + A) 2 
Hence 

q 1 / -iaiki-ia z k^-ivk-\ -f3m + a 2 k 2 \ , . 

A /c 2 +m 2 + (if /c + A) 2 \ (3m - a 2 k 2 -iaiki -ia 3 k 3 - ivk -\ J 

Taking the inverse Fourier transform we obtain 

„ ( v / («i9i + a 3 <9 3 + t> • V - A) -(/3m + a 2 <9 2 ) \ , x no r\ 

GM ={ (Pm + a 2 d 2 ) (a ldl + a 3 d 3 + v ■ V - A) ) ^ ^ 

where 

seR3 (13 - 7) 

Note that denominator in RHS (13.7) does not vanish for Re A > since \v\ < 1. This 
implies 

Lemma 13.1. The operator G\ with the integral kernel G\(y — y') , is continuous operator 
L 2 @L 2 ^L 2 @ L 2 for Re A > 0. 

From now on we use the system of coordinates in y-space in which v = (\v\, 0, 0), hence 
vk = \v\ki. Let us compute the function g\(y). One has 



k 2 +m 2 + (i\v\k 1 + \) 2 = ^k 2 + k% + kl + 2i\v\k 1 \ + \ 2 +m 2 = ^(k!+i-f 2 \v\\) 2 + kl + kl + x 2 
where 



7 = 1/VT^~ 2 , x 2 = -^- + A 2 +m 2 = -^— 2 +m 2 = 7 2 (A 2 +/i 2 ), p := m/ 7 (13.8) 

1 — v 1 — i> 

Hence formally, 

1 f e~ iky dk e -7 2 MA?,i r e - iky alk 



9\{y) 



(2tt) 3 / 2 / ^(jfei + i 7 2 |t;|A) 2 + fc! + fc| + k 2 " (2tt) 3 / 2 / ^2 A; 2 + k 2 + k 2 + n 2 
)e -i\v\\yi r e~ ik ydk 
(2tt) 3 / 2 ./ fcf + k 2 + k1 + k 2 



je / e dk = ie - llv \xy im _ k2) (13 0) 
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Here y\ = 73/1, y = (73/1,2/2,2/3), and i?(y — y',() is the integral kernel of the operator 
R(C) = (-A - Cr 1 - It is well known that R (y,() = / Ait\y\. Therefore, 

e - -x\y\->tiyi 

9x(y) = , I-, (13.10) 
47r|y| 

where xi := 7|f |A. We choose Reft > for Re A > 0. Let us note that for < \v | < 1 

0<Rex!<Rex, Re A > (13.11) 
Let us state the result which we have got above. 

Lemma 13.2. i) The function g\(y) decays exponentially in y for Re A > 0. 
ii) The formulas (13.10) and (13.8) imply that for every fixedy, the function g\(y) admits 
an analytic continuation in A to the Riemann surface of the algebraic function a/A 2 + /i 2 
with the branching points A = ±ip. 

Thus, from (13.1) and (13.2) we obtain the representation 

fx = -G 1 ^ Ql -G 1 ^ Q2 -{GfV Pl )-Q+{GfV P2 )-Q (13.12) 
^2 = -Gfm Ql -Gfm Q2 -{GfV Pl )-Q+{GfV P2 )-Q 

Step ii) Now we proceed to the last two equations (13.1): 

- \Q + B V P = -Q , (V^-, Q ■ Vp 3 ) - (W,, Pj ) - XP = -P . (13.13) 
Let us rewrite equations (13.12) as = ^ j(Q) + ^(^o), where 

*x(*o) = -G?* 01 - G\ 2 V 02 , tf 2 (* ) = -Gf^oi ~ Gf * 02 (13.14) 
MQ) = (-G\ 2 V Pl + Gl'Vpi) ■ Q, MQ) = (-GfV Pl + GfV P2 ) ■ Q (13.15) 
Then (V^j,pj) = (V^j(Q),pj) + (V^(* ), Pj), and the last equation (13.13) becomes 

(V^,Q • V Pj ) - (W,(Q),p,> -\P= -P + (V^(tf ),P;> =: -P - $(A) 

where 

*(\) = {* j (* ),Vp j ) (13.16) 

First we compute the term 

Applying the Fourier transform F y ^ k , we have by the Parseval identity and (A. 20) that 

^2{diip vj ,dipj) = y^^-ikj^j, -ikipj) = / hki{^ vl ■ p x + ip v2 ■ p 2 )dk (13.17) 
3 3 J 

- - [ k km ^ l ' ^ + ^ 2 ' p2 dk - - \ k i k i®( k ) dk 
J % k 2 + m 2 — (\v\ki) 2 J k 2 + m 2 — (\v\ki) 2 



--: -Li 



As the result, (Vip V j,Q ■ Vpj) = —LQ, where L is the 3x3 matrix with the matrix 
elements Lu. The matrix L is diagonal and positive defined by (1.10). 
Now let us compute the term —(V$?j(Q),Pj) = (^j(Q),Vpj). One has 

d iPj ) =J2 ((-G?d lPl +G^d lP2 , d iPl )-(Gfd lPl -Gfd lP2 , d iP2 ))Qi =J2 Hu(\)Qi 
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and by the Parseval identity and (1.3)-(1.4) we have 

Hu(X) : = {-G^d lPl + G 1 x 1 d lP2 ,d iPl )-{Gfd lPl -Gfd lP2 ,d iP2 ) 

= {[{(3m - a 2 k 2 )pi - (iaih + ia 3 k 3 + i\v\ki + X)p 2 ]g x k h k iPl ) 
+ {[{ia-tki + ia 3 k 3 + i\v\k! + X) Pl + {/3m - a 2 k 2 ) P2 \g x ki,k iP2 ) (13.18) 
gpi ■ Pi + Pp2 ■ 92 dk _ f hhB{k)dk 
k 2 + m 2 — {\v \ki — iX) 2 J k 2 + m 2 — {\v \ki — iX) 2 

The matrix H is well defined for Re A > since the denominator does not vanish. The 
matrix H is diagonal. Indeed, if i ^ I, then at least one of these indices is not equal to 
one, and the integrand in (13.17) is odd with respect to the corresponding variable. Thus, 
Hu = 0. As the result, (fyj(Q), Vpj) = HQ, where H is the matrix with matrix elements 
H a . Finally the equations (13.13) become 

»W ( p ) = ( P„ + %A) ) • M < A > = ( L -H(X) ~XE ) < 1319 ) 

Assume for a moment that the matrix M(A) is invertible (later we will prove this). Then 
we obtain 

(?) = "-'(*) (*,+%)• ReA>0 (1320) 

Finally, formula (13.20) and formulas (13.12), where Q is expressed from (13.20), give the 
expression of the resolvent -R(A) = (A — A)" 1 , ReA > 0. 

Lemma 13.3. The matrix function M(X) (respectively, M~ l (X)), ReA > admits an 
analytic (respectively meromorphic) continuation to the Riemann surface of the function 
vV + A 2 , A G C. 

Proof. The analytic continuation of M(A), exists by Lemma 13.1 ii) and the convolution 
expressions in (13.18) by (1.9). The inverse matrix is then meromorphic since it exists for 
large ReA. The latter follows from (13.19) since H(X) 0, ReA -> oo, by (13.18). □ 




14 Regularity on imaginary axis 

Let us describe the continuous spectrum of the operator A = A VjV on the imaginary axis. 
By definition, the continuous spectrum corresponds to w 6 1, such that the resolvent 
R(iu + 0) is not a bounded operator in S. By the formulas (13.12), this is the case 
when the Green function G\{y — y') loses the exponential decay. This is equivalent to the 
condition Re x = 0. Thus, ico belongs to the continuous spectrum if |a>| > p = my/l — v 2 . 
By Lemma 13.3, the limit matrix 

MM : =M(^ + 0)=( L _^ + 0) weR (14.1) 

exists, and its entries are continuous functions of ui G K, smooth for \oj\ < p and \oj\ > p. 
Recall that the point A = belongs to the discrete spectrum of the operator A by Lemma 
5.3, hence M{iu + 0) (probably) also is not invertible at u — 0. 

Proposition 14.1. (cf. [8, Proposition 15.1]) Let p satisfy the conditions (1.9)- (1.10), 
and \v\ < 1. Then the limit matrix M(iu + 0) is invertible for w / 0, w G K. 

Corollary 14.2. The matrix M^iiu) is smooth inu G M outside three points u> = 0, ±p. 
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15 Singular spectral points 

The components Q(t) and P(t) are given by the Fourier integral 

( )= — f e^M-^iu + O) ( D g ° . , ^ dco (15.1) 

if it converges in the sense of distributions. Corollary 14.2 alone is not sufficient for 
the proof of the convergence and decay of the integral. Namely, we need an additional 
information about behavior of the matrix M _1 (iu;) near its singular points 00 = 0, ±/i, 
and asymptotics at \uu\ — > 00. We will analyze all the points separately. 

I. First we consider the points ±/i. 
Lemma 15.1. The matrix M _1 (iu) admits the asymptotics in a vicinity of ±p: 



M-\iu) = C ± + 0((ut^), d u M~\iu) = 0{(cot^), ^M" 1 ^) = 0((u;^)-§) 

(15.2) 

Proof. It suffices to prove similar asymptotics for M(iu>). Then (15.2) holds also for 
M^^cu), since the matrices M(±ip) are invertible. The asymptotics for M(iu) holds by 
the convolution representation (13.18) 

Hjj(X) = (mgx/3 * djpi, djpi) + (mg x /3 * djp 2 , djp 2 ), j = 1, 2, 3 (15.3) 

since g x admits the corresponding asymptotics by the formula (13.10). Namely 

9\{y) = + r ±(\y), A — >■ iz^i, ReA>0 
47T\y\ 

where 

r ± {\y) = 0{{\Ti^), 9xr ± (A,y)=0((A=FVi)-*), S^r ± {\y)=0{{l+\y\){\^i^) 
The condition (1.9) provides the convergence of all integrals arising in d^Hjj. □ 

II. Second, we study the asymptotic behavior of M _1 (A) at infinity. 
Lemma 15.2. There exist a matrix D and a matrix-function Di(u>), such that 



M~ 1 (iuj) = — + D 1 (uj), M>/i + l, coeR (15.4) 
to 



where, for k = 0,1,2 



|8£DiM| < \u\ >fi + l, coeR (15.5) 



Proof. The structure (14.1) of the matrix M{iu) provides that it suffices to prove the 
following estimate for the elements of the matrix H(iu) := H(iu + 0): 

\d^H j:j (X)\<C(k), AeC, |A|>/i+l, j = 1,2,3, A; = 0,1,2 (15.6) 

The estimate (15.6) follows from the representation (15.3) and the bounds 



|<7a(s/)| < r-p |^A(y)|<^ + C 3 , \dlg x (y)\<-^ + C 5 \y\, ReA > 

1 2/ 1 I I 2/ I 



□ 
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III. Finally, we consider the point u = which is an isolated pole of a finite degree 
by Lemma 13.3. In Appendix B we prove that the matrix M -1 (?u;) can be written in the 
form 

M-\iu) = (15.7) 

where i, j — 1,2 are the diagonal matrices, smooth for the u G (— //). Moreover, 

M 22 = M U , M 11 = iM 12 B- 1 (15.8) 



16 Transversal decay for the linearized equation 

Here we prove Proposition (6.7). 

16.1 Decay of vector components 

First, we establish the decay (6.20) for the components Q(t) and P(t). 
Lemma 16.1. Let X G Z v fl £ v . Then Q{t), P(t) are continuous and 

\Q(t)\ + \p(t)\<c u ( P ,v)(i + \t\r 3 / 2 , t>o. 



16.1) 



Proof. The expansions (15.2), (15.4) and (15.7) imply the convergence of the Fourier 
integral (15.1) in the sense of distributions to a continuous function of t > 0. Let us prove 
(16.1). First let us note that the condition X G Z v implies that the whole trajectory 
X(t) lies in Z v . This follows from the invariance of the space Z v under the generator 
A v>v (cf. Remark 6.6). If X £ Z v , then the components Q(t) and P(t) may contain 
non-decaying terms which correspond to the singular point oj = since the linearized 
dynamics admits the secular solutions without decay, see (6.16). We will show that the 
symplectic orthogonality condition leads to (16.1). Let us split the Fourier integral (15.1) 
into three terms using the partition of unity Ci( w ) + (2^) + Ca( w ) = 1, w 6 K: 



( p(t) ) = 27 / e ^ (ClM + Uu) + CsH)^- 1 ^ + 0) 

where the functions Cj( w ) £ C°°(R) are supported by 



Qo 
P + 



3=1 



supp Ci c {to e 

supp ( 2 C {co G 

supp C 3 C G 
i) Let us represent Ij(t), j = 1,2 as 



e /2 < |w| < /I + 2} 
|u;| > ii + 1} 
M < e } 



16.2) 



1 
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Auit, 



Qiuj) M-\iu + 0) 





/(*) 



Qo 



+ M" 1 (iw + 0) 







du 



where 



Sj (t) = A-'^wjM-^iw + 0), /(t) = A _1 $(ia;) 



(16.3) 
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By (13.14) 

*i(*o) = -AReW+(t)(*io + i*2o), * 2 (*o) = -Mm W+(t)(V 10 + M 20 ) 

where W+(t) is the dynamical group of the equation 

V(x, t) = [a ■ V + i(3m + v ■ V]^(x, f) (16.4) 

Evidently, for the group the W+(t) the bound (17.5) obtained in Lemma 17.1 for the 
group W v (t) also holds. Hence, (13.16) implies 

|/(*)| = |Re(W+(0(*io + ^2o),Vp)| < C u (p,v)(l + t)- 3 / 2 (16.5) 

Further, the function Si(t) decays as (1 + |t|)~ 3 / 2 by asymptotics (15.2), and the function 
s 2 (t) decays as (l + |t|)~ 2 due to Proposition 15.2. Hence, formula (16.3) implies the decay 
(1 + |t|)- 3 / 2 for h(t) and I 2 (t). 

iii) Finally, the function I 3 (t) decays as if Zq E Z v . It follows from next lemma 
Lemma 16.2. If Z G Z v then 

(9M)=m-V)( Po + V)) eC "(-"-") (16 ' 6) 

Proof. In Appendix C we prove that the symplectic orthogonality conditions (6.7) at t — 
imply 

P + $(0) = 0, B^Qq + $'(0) = (16.7) 

Then 

P + = - $(0) = iwTiH 

where T 3 -(w) e C°°(-/x, //), since *j(* ) e C°°(-n,n), by (13.14) and (13.6). Therefore, 
representations (15.7)-(15.8) imply 

= M2i(w)Qo + «M22MTi(w)eC 00 (-/i,/ i ) 
Q(iw) = -M 11 (co)Q + -M 12 (co)T 1 (co) = -M 12 (B; 1 Q + T 1 (co)) 

U! CO LO 

= -M 12 T 2 (u)eC co (- f i,fi) 

□ 

16.2 Decay of fields 

Now we prove the decay of the field components ^i(x, t), ^ 2 (x, t) corresponding to (6.20). 
The first two equations of (6.15) may be written as one equation: 

V(x,t) = [-a-V-ipm + v-V]V(x,t)-iQ(t)-Vp, iGK 3 , teR (16.1) 

where ty(t) = ^i(-,t) + W 2 (-,t)). Applying the Duhamel representation, we obtain 

q?(t) = W v (t)V - f W v (t - s)Q(s) ■ Vp ds, t>0 (16.2) 
Jo 

where W v (t) the dynamical group (propagator) of the "modified" free Dirac equation 

y(x, t) = [-a ■ V - ifim + v ■ V]*(x, t) (16.3) 

Lemma 17.1 on the weighted decay for the group W v (t), the decay of Q from (16.1), and 
representation (16.2) yield 

11^(011- < C v (p,v)\\V \\ v (l + |t|)- 3/2 , * > (16.4) 
for any ty & Z V C\ E v . It completes the proof of Proposition 6.7. 
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17 Weighted decay for free Dirac equation 

Lemma 17.1. For any $6Lj with v > 3/2 the bound holds 

\\W v (t)n-„<^^, t>0 (17.5) 

Proof. Step i) Note, that 

(d t + a ■ V + i(3m + v ■ V) (d t - a ■ V - i(3m + v ■ V) = (<9 2 + 2<9 t v • V - A(l - w 2 ) + m 2 ) 

Hence the integral kernel W v (x — y,t) of the operator W v (t) has the form 

W v (z, t) = {d t + a-V + i(3m + v-V) G v (z, t), (17.6) 

where G v (z,t) is a fundamental solution of the "modified" Klein-Gordon operator 

(<9 2 + 2d t v • V - A(l - v 2 ) + m 2 )G,(2, t) = S(z)S(t) 

Let G v (t), t > be the operator with the integral kernel G v (x — y,t). It is easy to see 
that 

[G v (t)Q](x) = [G (t)$\(x-vt), xeR 3 , t>0 

Then 



G.(,,t) = <*,(,-*, t) = j(< "J' = " '* ~ "pt^ ~ '* - Jff>, ( >0 

47Tt 47T ^2 — \z — Vt\ 2 

where J\ is the Bessel function of order 1, and 9 is the Heavyside function. Let us fix an 
arbitrary e e (|t>|, 1). Well known asymptotics of the Bessel function imply that 

\d t G v {z,t)\, \d Zj G v (z,t)\ <C(e)(l + t)- 3 / 2 , \z-vt\<et, t>l, j = 1,2,3 (17.7) 

Step ii) Now we consider an arbitrary t > 1. Denote = e — \v\. We split the function 
<3> in two terms, $ = $i ;t + $ 2 ,* such that 

||$MllLg + ||$2, t |Ug <C||$|| L 2, *>1 (17.8) 

and 

$ 1>t (x) = for \x\ > and $2,t( x ) = for M < ~7~ ( 17 - 9 ) 

The estimate (16.4) for W v (t)$2,t follows by charge conservation for Dirac equation, (17.9) 
and (17.8): 

„ I. . . | C(e)\\$2t\\v C' 1 (£)||$|| L 2 

(17.10) 

since z/ > 3/2. 

S'tep raj Next we consider W v (t)$i jt - Now we split the operator W v (t) in two terms: 

w v (t) = (i - c)w w (t) + cw„(f), * > i 
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where ( is the operator of multiplication by the function such that ( = ((s) e 

Cg°(R), C(s) = 1 for |s| < ex/4, £(s) = for \s\ > e l /2. Since 1 - C(MA) = for 
\x\ < Sit /A, then applying the charge conservation and (17.8), we have for t > 1 

ii(i-o^(o^ii^< — — = {1+t)u < {1+t)u < (1+t)3/2 

(17.11) 

since z/ > 3/2. 

iSfep w) It remains to estimate (W v (t)$ij. Let x* be the characteristic function of the 
ball |x| < Sit/2. We will use the same notation for the operator of multiplication by this 
characteristic function. By (17.9), we have 

CW v {t)^ t = (W v (t) X t® (17.12) 

The matrix kernel of the operator (W v (t)xt is equal to 

W' v {x -y,t) = ((\x\/t)W v (x - y, t) Xt {y) 

Since C(MA) = for \x\ > s x t/2 and Xt(y) = for \y\ > s x tj2. Therefore, W' v (x-y,t) = 
for \x — y\ > sit. On the other hand, \x — y\ < s\t implies \x — y — vt\ < st, since £i + |i>| = e 
by definition of s±. Hence, equality (17.6) and bounds (17.7) yield 

\W^x-y,t)\<C(l+ty 3 / 2 , t>\ (17.13) 

The norm of the operator (W v (t)xt '■ L 2 — > L 2 _ v is equivalent to the norm of the operator 

(x)- v CW v (t)xt(y)(y)- v :Ll^Ll (17.14) 

Therefore, (17.13) implies that operator (17.14) is Hilbert-Schmidt operator since v > 3/2, 
and its Hilbert-Schmidt norm does not exceed C(l + 1)~ 3 / 2 . Hence, by (17.12) and (17.8) 

\\(W v (t)<5> ht \\ L2 _ v < C(l+t)- 3 / 2 \m Lh t>l (17.15) 
Finally, the estimates (17.15), (17.11) and (17.10) imply (16.4). □ 

A Computing Q(ri, Tj) 

Here we justify the formulas (3.7)-(3.9) for the matrix Q. 

1) First, the Parseval identity implies 

Q(Tj, ti) = (djip vl , dii> v2 ) -{djip v2 , diip vl ) = J kjki dk(^p vl -^ v2 -^ v2 -^ v i) = 0, j, I = 1, 2, 3 

since the integrand is odd function. 

2) Second, we consider 

fi(r j+3 , r l+3 ) = (d Vj ip vl , d Vl ip v2 ) - {d Vj ip v2 , d Vl ip vl ) (A.16) 
Let us derive the formulas for ip v i and tp v2 . The first equation of (2.4) implies 

[{vjdj) 2 — A + m 2 ]ip v = [ivjdj + iajdj — j3m] p\ 

Hence 

[( v jdj) 2 - A + m 2 }i[) v i = -[vjdj + aidi + a 3 d 3 }p 2 - [a 2 d 2 + fim\p x 
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[(vjdj) 2 - A + m 2 \i\) v2 = [vjdj + + a 3 d 3 }pi - [a 2 d 2 + /3m]p 2 
Applying the Fourier transform, we obtain 

[ivk + iaiki + ia 3 k 3 ]fa + [a 2 k 2 - f3m]fa 



1pv2 = 

Differentiating, we get 

d Vj 4> v i = 

Hence, (A. 16) implies 



— (v k) 2 + k 2 + m 2 
■[ivk + iaiki + ia 3 k 3 ]fa + [a 2 k 2 - /3m]p 2 

— (vk) 2 + k 2 + m 2 



(A.17) 



ikjp 2 



2kjvkip v i 



-(vk) 2 + k 2 + m 2 
-ikipi 



+ 



-(vk) 2 + k 2 + m 2 
2k x vk^ v2 



vk) 2 + k 2 + m 2 —(vk) 2 + k 2 + m 2 



(A.18) 



,'kjkilp! ■ fa- fa- fa]dk f Akjki(vk) 2 [tp v i ■ tp v2 - i) v2 ■ ip vl ]dk 
M{Tj+ 3 , n +3 ) = I — , lxoxo v 



+ 



I 



(k 2 + m 2 — (vk) 2 ) 2 J (k 2 + m 2 — (vk) 2 ) 2 

2ikjkvk[p 2 ■ V>„ 2 + 4>v2 ■ fa + fa • i>vi + 4>vi • fa] 



(k 2 + m 2 — (vk) 2 ) 2 







since all integrands are odd functions. 
3) Finally, (A.18) implies 

Q(rj, t 1+3 ) = -{djip vl , d Vl i) v2 ) + {djiprt, d Vl tp vl ) + ej ■ d Vl p v 

ikjtp vl ■ [-ikipi + 2kivkip v2 )] - ikji) v2 ■ [ihfa + 2kivkip vl )] 



kjki 



k 2 + m 2 — (vk) 2 
-$vi ■ fa + 4>v2 ■ fa] + 2ivk[tp vl ■ ip v2 - i) v2 ■ $ vl ] 



(A.19) 

dk ~\~ 6j ■ du^Pi) 



dk ~\~ Gj ' du^Pi) 



k 2 + m 2 — (vk) 2 

Recall, that Pj(x) are even, then fa(k) are real. Hence (1.3)-(1.4) and (A.17) imply 

(k 2 + m 2 - (vk) 2 )(i[j vl ■ fa + ip v2 ■ fa) = [a 2 k 2 - (3m]fa ■ fa + [a 2 k 2 - f3m]fa ■ p (A.20) 
+ [ivk + iaiki + ia 3 k 3 ]fa ■ fa — [ivk + iaiki + ia 3 k 3 ]fa ■ fa = —Bp ■ p 



(k 2 + m 2 - (vk) 2 ) 2 (ip v i ■ 4> v2 - ^ v2 ■ $ vl ) = 2i(k 2 + m 2 - (vk) 2 ) 2 lm$ vl ■ 4> v2 ) (A.21) 
= —2f3mfa ■ [ivk + iaiki + ia 3 k 3 ]fa — 2[ivk + ia\ki + ia 3 k 3 ]fa ■ f3mfa = —2ivkBp ■ p 

Substituting (A.20) and (A.21) into the right hand site of (A.19), we obtain 

B(k) A(vk) 2 B(k) 



tt(Tj,T l+3 ) = J kjki^ 



+ 



.(k 2 + m 2 - (vk) 2 ) 2 (k 2 + m 2 - (vk) 2 ) 3 
that correspond to (3.7) - (3.9). 



^jdk + Cj ■ d Vl p v 
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B Computing M (iuj) 

Here we derive formula (15.7). Denote F(u) := —L + H{iu + 0) which is diagonal. Then 
by (14.1) for u G K we obtain 



detM(ioo) = det 
where 



iuE 
-F(u>) 



iuE ) (' 



u 2 + 



F 22 (oo) 



7 



)(- 2 +^) (B.22) 



k 2 Bdk ( - - 

\ m 2 + k 2 — (\v\ki + iu) 2 m 2 + k 2 — (\v\ki) 2 



, j = 1,2,3 (B.23) 



Formula (B.22) is obvious since both matrices F(u) and B v are diagonal, hence the matrix 
M(iu) is equivalent to three independent matrices 2x2. Namely, let us transpose the 
columns and rows of the matrix M(iu) in the order (142536). Then we get the matrix 
with three 2x2 blocks on the main diagonal. Therefore, the determinant of M(iu) is 
product of the determinants of the three matrices. Further, 

M 21 (u) M 22 (uo) 



where 



M u (ou) = M 22 (u) 



M-\iuj) = 
( 

V 



(B.24) 



w 2 7 3 +Fn(w) 










w 2 7+F22(w) 










\ 



M- 



12 












-1 











-1 

cj 2 7+F 33 



/ -7^11 
' w 2 7 3 +*ii 



M- 



21 



V 










W 2 7+F22 










\ 



- 7^33 i 
J 2 7+^33 / 



Let us prove that for u G (— /x, /x) 

F jj ( U )=u 2 f jj (u), /^HgC^-^/x), f jS (0)>0 (B.25) 
Indeed, formula (B.23) implies that Fjj(0) = 0. Differentiating (B.23), we obtain 



k 2 B{k)dk 







fci 


[k 2 + m 2 - 


-(\v 


h) 2 ) 2 



since integrand is odd function in respect to k±, and 

k 2 + m 2 + 3(\v\ki) 2 



F;;.(0) = 2 J k 2 B{k)dk-^ 



(k 2 + m 2 - (\v\hy) 3 
By (B.25) we can represent the matrices JWjj(w) as 



= 



> 



Mn(w) = M 22 (w) = - 



V 



-«7" 



7 3 +/ll(^) 








— vy 



7+/22M 










\ 



7+/33(w) / 



M 12 (w) = 



V 



7 3 +/nH 










-1 



/ -7 3 /nM 
' 7 3 +/nH 



M 2 i(w) = 
where A^ij(cu) G C°°(-/x, /x) 



V 



o 
o 



7+/22M 




-7/22 M 
7+/22M 










\ 



7+/33M / 




(B.26) 



\ 





-7/33 M 

7+/33(w) 



-M 2 i(cu 
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C Symplectic orthogonality conditions 

Here we derive conditions (16.7) from the symplectic orthogonality conditions (6.7). First 
let us compute $(0). Formulas (13.14) and (13.16) imply 

($(0)), = (G^oi + G 12 i> 02 ,ik j p 1 ) + (G^ 02 - G l2 ^ m ,ik 3 p 2 ), j = 1,2,3 

On the other hand, by (13.5) formulas (A. 17) read 

= -G^h + G 12 Pl , ^v2 = Gl 1 ^ + G l2 p 2 

Hence, for j = 1,2,3 

= -n(Z Q , Tj) = (* i, dtfa) - (*02, djil> vl ) + P • e 3 

-(^WM^oVi + G l2 p 2 )) + ($02, ikj(Gl 2 pi - Gl l p 2 )) + P • e, = ($(0) + P ) j 

since (G* 1 )* = -Gj 1 , (G* 2 )* = G* 2 . Hence the first condition (16.7) follows. Further, 



^Gi 1 



-1 - 2ivkG$ 



-2ivkGl 2 



a=o A; 2 + m 2 — (vk) 2 ' a=o k 2 + m 2 — (vk) 2 

Then (13.14) and (13.16) imply for j = 1,2,3 

^mh / j> m + 2ivk{G^i> m + G 12 i> 02 ) A \ / ^^^(G^-G^qx) A 
($ (0)) ^-\ A: 2 + m 2 -(,A:) 2 '^r\ fc 2 + m 2_ W 

On the other hand, from (A. 17) and (A. 18) it follows that for j = 1, 2, 3 

ikjp 2 + 2k j vk(-G 1 ^p 2 + Gj 2 pi) ? _ -ifcj-pi + 2/c i wA;(Gj 1 pi + Gj 2 p 2 ) 



" J " A; 2 + m 2 — (t>A;) 2 " J " A; 2 + m 2 — (t>A;) 2 

Hence, 

= Q(Z , r j+3 ) = (^oi, ^.^ 2 ) - (*02, d Vj ip vl ) + Q • 9 v .p v 

. T , -ik jPl + 2fc j t;/c(G^ 1 p 1 + G* 2 p 2 ) \ / ik jP2 + 2k j vk(-G 1 1 p 2 + G 12 ^) 



k 2 + m 2 — (vk) 2 I \ ' A; 2 + m 2 — (t>A;) 2 

+ Qo • ^.p, = ($'(0) + B^Qoh, j = 1,2,3 

since Q ■ d Vj p v = Q ■ B~ x e^ = B~ x Qq ■ ej. Hence the second condition (16.7) follows. 
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